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ABSTRACT 

An analysis of the nonlinear development of the large-scale 
structures or instability waves in compressible round jets 
was conducted using the integral energy method. The equa- 
tions of motion were decomposed into two sets of equa- 
tions; one set governing the mean flow motion and the 
other set governing the large-scale structure motion. The 
equations in each set were then combined to derive kinetic 
energy equations that were integrated in the radial direction 
across the jet after the boundary-layer approximations were 
applied. Following the application of further assumptions 
regarding the radial shape of the mean flow and the large 
structures, equations were derived that govern the nonlin- 
ear, streamwise development of the large structures. Using 
numerically generated mean flows, calculations show the 
energy exchanges and the effects of the initial amplitude 
on the coherent structure development in the jet. 

INTRODUCTION 

Strong evidence suggests that jet noise, particularly in 
the supersonic regime, is dominated by contributions from 
wave-like structures in the initial region in the jet that are 
associated with Mach-wave radiation. These structures 
cannot be captured by classical turbulence modeling. Di- 
rect Numerical Simulations (DNS) and Large Eddy Simu- 
lations (LES) can successfully capture these structures, but 
they are computationally intensive. 

There have been several attempts to model this structure 
using the integral energy approach. The basic idea is that 
the coherent structure is modeled as nonlinear instability 
waves that grow and decay along the jet. Each flow com- 
ponent is split into three quantities; a time-average mean 
component, a large-scale wave-like coherent component, 
and a fine-scale random turbulence component. Starting 
from the full compressible Navier-Stokes equations and us- 
ing time averaging, the momentum equations for each flow 
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component can be obtained. These equations are then used 
to derive the kinetic energy equations, which are integrated 
across the jet to produce a set of simultaneous, ordinary 
differential equations describing the interactions among the 
various scales of motion. Shape assumptions are then made 
to close these equations to obtain the development of each 
flow component. The plane shear layer case has been stud- 
ied by several researchers. 1,2 The incompressible round 
jet case has been studied for single and multi-frequency 
modes. 3,4 The compressible planar shear layer case has 
been studied by Lee & Liu. 5 The focus of this study is 
the compressible, supersonic round jet case because of its 
obvious relevance to practical applications. 

Apart from modeling the physics involved, the study 
of the development of the coherent mode was conducted 
to simulate excited jets. Excitation via single or multi- 
frequency mode has been conducted in the past to control 
the jet signature. For instance, Raman et al. 6 show that by 
exciting an incompressible round jet via forced waves of 
various frequencies, amplitudes, and phases, the jet can be 
made to spread faster and thus may reduce its noise signa- 
ture. Arbey & Ffowcs-Williams 7 considered the objective 
of reducing the peak noise by forcing an excitation wave at 
the subharmonic frequency of the peak noise. They showed 
that an interaction process could occur that would result 
in suppressing the fundamental and thus reducing the peak 
noise. 

The starting point of this analysis is that we consider 
a turbulent round jet at a sufficiently high speed so that 
the compressibility is significant. The development of this 
jet in the unexcited case is assumed to be given by some 
other means (e.g. analytically or via Reynolds averaged 
numerical simulations). This jet is then excited by a single- 
frequency instability wave. The nonlinear development of 
this wave will be presented herein based on the integral 
energy approach. Along with the wave development, the 
mean flow-spreading rate is also modified and will be pre- 
sented herein. Since the focus here is the supersonic jet, we 
must consider the helical modes as they are more amplified 
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than the axisymmetric ones. 8 Once the computations are 
completed, a description of all the variables of the large- 
structure disturbance field in and near the jet is given for 
a particular mode number and frequency. It can be used 
as the acoustic source to compute the noise radiation field 
outside of the jet. This is described in a separate paper. 9 

FORMULATION OF THE PROBLEM 

Consider a high-Reynolds number turbulent jet issuing 
from a nozzle of diameter D in a still air. The jet is 
shock-free, but the Mach number is high enough for com- 
pressibility effects to be significant. The jet is excited 
by a single-frequency instability wave of Strouhal num- 
ber fD/Uj. The density and the component velocities are 
normalized by the jet exit density and axial velocity at the 
centerline, pj and Uj, respectively. The pressure is normal- 
ized by PjUj , time by D /U j, and spatial coordinates by 
D. Each flow parameter is split into a time-averaged part 
Ui(x, r, (f>), a coherent, wave -like part, u'^x, r, (f>, t), and a 
fine-scale random turbulence part, u"(x, r, <j>, t). Thus, the 
velocity can be written : 

Ui = Ui(x,r,<j>) + u' i (x,r,<j>,t) + u”(x, r, <j), t) (1) 

where i = 1,2,3. In the cylindrical coordinates, 1 refers 
to the axial direction x with axial velocity u, 2 refers to the 
radial direction r with radial velocity v, and 3 refers to the 
azimuthal direction (j> with azimuthal velocity w. An over 
bar, ( ), denotes a time-averaged quantity. The pressure 

and the density are similarly split: 

p = P{x, r, <t>) + p'( x, r, cj), t) + p"(x, r, cj), t) (2) 

P = P(x, r, cj)) + p'(x, r, cj), t) + p"(x, r, cj), t) (3) 

In the subsequent analysis, however, the random tur- 
bulence component will not be explicitly considered any 
further, thus we have 

u i = U i (x,r,<j)) + u' i (x,r,^,t) (4) 

P = P(x,r,cj>) + p'(x,r, cp,t) (5) 

p = p(x,r,4>)+p'(x,r,<j),t) (6) 

For the product of the density with the velocities, we get 

(pui) = (p + p')(U i + u'j) = pU i + pu' i + p'U i + p'u'i 
Time-averaging yields 

(Wi) =pUi + p'u'i ( 7 ) 

and we define 

Ui = piLi — pui =pu'i+p l U i + p'u'i - p'u'i . 

Neglecting the difference between the product of distur- 
bance variables and the averaged value of the product re- 
sults in 

Ui « pu'i + p' Ui. (8) 


Equations of motion 

The formulation begins with the following nondimen- 
sionalized continuity and momentum equations in cylindri- 
cal coordinates: 

Pt + ( pu) x + ^ [r(pv)\ r + ^(pui)# = 0 (9) 


(pu) t + (pu 2 +p) a 


( 10 ) 


+ \ V (puv)\ r + ^ ( puw ) 4> = ^- A u 


(. pv)t + (puv) x (11) 

+ 1 [r (p + pv 2 )] r + 1 (. pvw) (j> - — 


P + - 
r Re 


Av j ( v + 2 w <t > ) 


(pw) t + ( puw ) a 


( 12 ) 


If/ M 1 / 2 \ puw 

+ -[r(pvw)] r +-(p + pw ) 4> + — 


1 

Re 


Aw - (w - 2v ( f > ) 


where the Reynolds number is Re = pjUjD/p. The sub- 
scripts denote differentiation and the Laplacian is 

d 2 1 d d Id 2 

dx 2 ^ r dr dr + r 2 d(f> 2 

Substituting relations (4) to (8) into (9) to (12), we get: 


(P + P')t + ( P u + u) a 


(13) 


+ ^ [r (pv + u)] r + ^ (pw + w)^ — 0 


(pu + u) t + [P + p' + (U + u') (pu + «)] x 
+ - \r (V + v') (pu + u)] r + ^ [ (W + w') (pu + u)] 

= T a(E7+ „-) 


(14) 


(pv + v) t + [(U + u') (pv + v)\ x 

+ ~[ r ( P + P' + ( V + V ') (PV + ^))] r 
+ - \(W + w') (pv + ti)] ^ — - (W + w ') (pw + w) 


1 


= -(P + P 0 (15) 

\^(V + v') -±(V + v> + 2(W + w') <p ) 


r 

1 

Re 
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(pw + w) t + [(U + v!) (pw + w)] x 
H — [r (V + v') (pw + w)] r H — (V + v') (pw + w) 
+ -[P + p' + (W + w')(pu + u)\ (16 


= Se [ A ( W + W -^( W + u/ ~ 2 ( V + 

These equations will be used to derive a set of equations 
governing the mean flow and a set of equations governing 
the large-scale structures. 

Mean flow 

The mean flow equations are obtained by time-averaging 
equations (13) to (16) containing the two-component de- 
composition. The continuity equation for the mean flow 
is 

m x +\[rm\r+l(m+= ° ^ 

and the three momentum equations for the mean flow are: 


(y U 2 + u'u + P 


) +1 K 

/ x r L V 


pUV + v' 


+ -(pWU + w'u) = — AU (18) 
r V / <j> Re 


(yUV + u’v'j H — r (yV 2 + v'v + P^j 
+ - (y W V + w'v'j — - (y w 2 


H — (y W V + w'v'j (y W 2 + w'w'j 

= t p+ s[ aF “^< F+2W *)] <19) 


(yUW + u'w'j H — r (y V W + v'w 


1 /_ f -'Z — — — 

-t — ( p W + w'w + P 
r \ 


) + - ( pWV + v'w 
/ s r \ 


r \ / <f> r \ / 

=1[ aw -M w - 2V ^\- (20) 

Large-scale structure 

The equations governing the large-scale structures are 
obtained by subtracting the mean flow equations (17) to 
(20) from equations (13) to (16). The second order terms 
are neglected with respect to first order terms and the terms 
u'jUj — u'plj are also neglected. The continuity and mo- 
mentum equations for the large-scale structures are: 

Pt + u x + - (rv) + -w^ = 0 (21) 

r r 


u t + (p 1 + pUu' + u U) x + - [r ( pUv ' + u V)] ] 
+ \ (pUw 1 +flW)^= ^ Au' 


v t + ( pVu ' + v U) x + - [r ( p ' + pVv' + v V)] ? 
H — (pVw 1 + vW)^ — - (jiWw' +wW) 

= V + ^ Av ' - y( v ' + 2w '<p) 


w t + [pWu 1 + w U) x + - [r ( pWv 1 + w V)] r 
4 — (p' + pWw' + w W) , + - (pWv 1 + w V) 


= is [ A ”' -M w '- 2v *( 


Kinetic energy equations 

Mean flow kinetic energy equation 
The mean flow kinetic energy equation is obtained by 
first multiplying the x-momentum equation (18) by U, the 
r- momentum equation (19) by V, and the ^-momentum 
equation (20) by W. Then, the resulting equations are 
added together. This combined equation is composed of 
several main terms including the stress terms, the pressure 
terms, and the viscosity terms. Using the mean flow con- 
tinuity equation and defining the mean flow kinetic energy 
2 2 2 

as K = (U + V + W )/2, the combined equation is re- 
arranged to obtain the kinetic energy equation for the mean 
flow 


— (y U K + u'u U + u'v V + u'w W + U P^J 

4 — — |V (y VK + v'u U + v'v V + v'w W + V P^j j 

4 — (pWK + w'u, U + w'v V + w'w W + W p\ 

r d(p \ ) 

-r(v,+ \(rV) r + \w^ 

— U x u'u — V x u'v — W x u'w — U r v'u 

— V r v'v — W r v'w U d>w'u V a,w'v 

r r 

1 ^ 1 ^ 1 

W <hw'w V w'w H — W v'w 

ak - 


10 / — 


-yA v + 2W «)-™(w-2vy 


where (U i X ) 2 = ( U x ) 2 + ( V x ) 2 + (W x ) 2 and similarly for 
(U ir ) 2 and (Uitf,) 2 . 

Large-scale kinetic energy equation 
The large-scale energy equation is obtained by multiply- 
ing the x-momentum equation (22) by u', the ?’-momentum 
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equation (23) by v' , and the ^-momentum equation (24) 
by w'. Again, the resulting equations are added together 
and the time-average is applied. This equation requires 
extensive manipulation to obtain its final form. The large- 
scale continuity equation (21) and equation (8) are used to 
combine terms. Defining the large-scale coherent structure 
kinetic energy as Q = ( u ,2 + v' 2 + w ' 2 ) /2 and noting that 
the mean flow is steady, the resulting energy equation for 
the large-scale structures is 

A®r 5 ) + i| (r? VQ) + \^WQ) 

W 'p'ch 

+ u'p' + v'p' H b u'u U x + u'v U r 

r 

u'w — — ~ — — — — v'w — — — — 

H U $ + v’u V x + v'v V r H V ^ + w'uW x 


+ w'v W r H Wa, H V H W 

r r r 


1 _ (yf )2 

= - AQ — (u' ix ) 2 — («' r ) 2 — — 


V w 


INTEGRAL FORM OF THE ENERGY EQUATION 

For the round jet, the mean quantities are assumed to 
be axisymmetric. Thus, W = 0 and d( )/d(j> = 0. The 
energy equations (25) and (26) for the mean flow and for 
the large-scale structures, respectively, become: 

(pUK + VulJ + v/vV + Up) 

H — r (~pVK + v'u U + v'v V + V p] 
r or L V / J 

-p(c7 s + I(rF) r ) (27) 

— U x u'u — V x u'v — U r v'u — V r v'v — -V w'w 

r 

= le\ AK ~( U ^ 2 -^r) 2 -^ 


d 10 vj’p’ 

-w~(pU Q) + - — (rpV Q) + u'p' x + v'p' r H 

ox r or r 

+ u'uU x + U ,: V U r + v'uV x + v'vV r + ^~^~V 

r 

i r _ ( u ' )2 

= - aq-( u 'j 2 -«) 2 -^F- c 


V w 

~ r 2(v' + 2 «^) - ^2 K - 2 ^) 


The usual boundary layer-type approximations are ap- 
plied to the mean quantities. The radial velocity is 
much less than the axial velocity, V ■C U, and the 
axial gradients are much less than the radial gradients, 
d( )/dx d( ) /dr. Second order terms are also ne- 
glected with respect to first-order terms. Thus, the mean 
flow energy equation reduces to: 

^^UPup) + ll[r(\- p VU^ + Vp)[ 

- P (u x + ^(rV)Jj -U r v'u (29) 


- ( u r y 


Re r " 


Multiplying equation (29) by r and integrating over r, we 
get: 

IS: 

/* OO POO 

— / PU x rdr— / P(rV) r dr 

Jo Jo 

/•OO y POO 

— / U r v'urdr = — — — / (U r ) 2 rdr 

Jo R-e J o 


2 -i r — »oo r / 2 \ ~\ 1 — »oo 

rpV — = 0 and r — 1 = 0 

2 \ 2 

J r —0 L \ / r J r —0 


The mean flow pressure is assumed to be constant across 
the jet. The integral form of the mean energy equation re- 
duces to the following simple form. 


^ poo ^ ^ r oo y r oo 

— -pU rdr = / U rv'urdr— — / ( U r ) 2 rdr 

dx J o 2 J 0 Re J 0 

(30) 

Similarly, the energy equation for the coherent compo- 
nent reduces to: 

1<PUQ )+ I^VQ) 

w'p\ 

+ u'p' + v'p' H + u'v U r 

r 

1 ( yj \ 2 

= Bi A <3 - KJ 2 - (O’ - -Jr- 


-J2( V ' + 2W '^-J2( W '- 2 V S) ( 31 > 


Multiplying (31) by r and integrating over r, we get: 


— ( pU Q)rdr = 
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into A(x). The radial shape functions are normalized such 
that 

(|u| 2 + |6| 2 + \w\ 2 ) rdr = 1 

Substituting the mode definitions into equation (30) and 
extending the definition of the averaging process to include 
the azimuthal direction. 



since 

[rpvm r o °°=0 and [rQ P ] ™ = 0. 

The physical interpretation of the terms appearing in the 
energy equations is obvious. For the mean flow equation, 
the left side is the mean flow advection of the mean flow 
kinetic energy. The first term on the right-hand side of (30) 
is the energy transfer from the mean flow to the coherent 
structure, and the second term is the viscous dissipation of 
the mean flow energy. As for the coherent mode equation 
(32); the left-hand side is the mean flow advection of the 
coherent mode kinetic energy. The first term in the right- 
hand side is the work done by the coherent mode pressure 
gradients, the second term is the energy transfer from the 
mean flow to the coherent mode, and the last term is the 
coherent mode energy dissipation. 

SHAPE ASSUMPTIONS 

To solve the above system of two energy equations, 
shape assumptions need to be made regarding the radial 
profiles. The mean flow is assumed to be a known function 
of r and 9, where 9 is the momentum thickness of the jet 
shear layer, and the coherent structure profiles are assumed 
to follow that of the locally-parallel, linear stability theory. 

Single azimuthal modes 

Considering single-frequency, single-azimuthal number 
modes, the normal mode composition suggests that these 
modes can be represented as: 

u'i(x, r, </>, t) = Ui(r)A(x) exp (*F + imp) + cc (33) 
p'(x, r, cp, t) = p(r)A(x) exp (vF + incp) + cc (34) 
p'(x, r, </>, t) = p{r) A(x) exp (vF + in<p) + cc (35) 

r 

*F = i / ad\ — iut 

where (") denotes the radial shape function of the trans- 
verse coordinate r at a given location along the jet. Ui(r), 
p(r), and p(r) are eigenfunctions corresponding to a given 
n and u>. Here, n is the azimuthal wave number indicating 
the rotation around the jet centerline, a is the wave number 
and cc denotes the complex conjugate. A(x) is the complex 
amplitude function of x and is to be determined by a non- 
linear analysis. In the nonlinear analysis, the linear growth 
rate as determined by the imaginary part of a is absorbed 


( ) = lim 

V ' T— >oo 27 tT 


( ) dt d<f>, 


where T is a large time scale, we obtain 

~7 ~T7T = \A\ 2 (I mw l + Imw 2 ) — 75 Imd, (36) 

dx dO Re 


where 


lam =2 pU r dr 


_dU 

I mw 1 = / (vu* + cc)p— rdr 

Jo 9r 

Imw2 = J (vp* + cc)U^rdr 


Im.rl. 


r°° fdu \ 2 


* rdr 


In this analysis, the mean flow is characterized by the mo- 
mentum thickness rather than by the axial distance. Hence, 
the integrals in the energy equations are dependent on 9 
instead of x. 

Upon substitution into equation (32), we obtain for the 
coherent mode: 


d{l aw \A\ 2 ) 


— 1^1 ( Imwl Imw 3 Ip Iwd 


/»oo 

I aw = (\u\ 2 + \v\ 2 + \w\ 2 )pUrdr 
Jo 

rv,* ^du , 

Imw 3 = J (Up + cc)V —r dr 


Iind — 2 


^/. W f. A\ 5k 7 

u(iap) 1 ( mp ) + cc rdr 

or r 


+ M(H 2 + H 2 + H 2 ) 


3 — 2 (l^l 2 + l^l 2 ) ^(’i-nvw* + cc) 


i~v Z r\ a / rv a Z \ 

ou ov ow 

TT + a- + rdr 

or or or 
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Pair of oblique waves 

Experimental observations have indicated that the helical 
modes appear in the form of a symmetrical pair with n, the 
mode number, set to both ±1 . 8 In this case, the normal 
mode definitions are modified as follows: 


u'( x, r, <t>, t) = u(r)A{x) exp (\k) cos (n<j>) + cc (38) 

v'(x, r, (/), t) = v(r)A(x) exp (\P) cos (n</>) + cc (39) 

w'(x, r, <f>, t ) = w(r)A(x) exp (\P) sin (n<j>) + cc (40) 

p'(x, r, </>, t) = p(r)A(x ) exp (\P) cos (n<j>) + cc (41) 

p'(x, r, <j>, t ) = p(r)A( x) exp (\D) cos(rn^) + cc (42) 

When these definitions are substituted into the integral 
equations (30) and (32) and the averaging process is ap- 
plied, we obtain the following sets of equations. 

~ dI ^ = \A\ 2 (Imwl + Im W 2) ~ Imd (43) 
ax dO Re 

where 


1 3 

lam =7, pU r dr 


f°° fdU\ 2 


Imd= L rdr 


d(l aw \A\ 2 ) 

dx 

where 





The radial shape functions are normalized such that 

l (N 2 + N 2 + N 2 )^r={^ = ° 

For the incompressible case, the above equations reduce to 
that of Lee & Liu . 10 


Mean flow 

In previous studies using the energy integral technique, 
the shape assumption for the mean flow was provided by 
analytic equations describing the velocity and density pro- 
files . 1,310 This works reasonably well in the initial po- 
tential core region of the jet where the shear layer is well 
described by a hyperbolic tangent or other functions and 
the locally-parallel flow instability wave is predominantly 
growing. In order to describe the full axial development 
of growth and decay of the instability wave, the mean 
flow must not only be described through the potential core 
region, but also possibly downstream into the fully devel- 
oped flow region. An analytic function description of the 
mean flow would not be smoothly continuous throughout 
this domain. Furthermore, an analytic function description 
of the mean axial velocity (e. g. hyperbolic tangent) com- 
bined with the Crocco-Busemann description of the density 
profile provides a nonphysical radial velocity profile. A 
proper physical description for the radial velocity profile 
is required for computing integral I mw 3 . These objections 
are overcome by determining the mean flow numerically. 
We compute the mean flow using the procedure of Dahl & 
Morris 11 with the addition that the radial velocity profile is 
computed from the mean continuity equation (17). Thus, 
the numerical calculations provide the mean flow profile 
shapes that smoothly transition from the initial region to 
the fully developed region and that can be parameterized 
by 9 as required by the analysis. 

The computational scheme for the mean flow of a jet was 
developed using a set of compressible, Reynolds -averaged, 
boundary layer equations with a modified mixing length 
model to determine the Reynolds stresses. In nondimen- 
sional form, the Reynolds stresses were computed by 


— pu'v ' = 


Pt dU 
PjUjD dr 


(45) 


where /it is a turbulent viscosity and is a function of r and 
x. Here, the unsteady terms represent all scales of turbu- 
lence. If the fluid viscosity is given by p, then in the shear 
layer /i T » p. 


Linear stability 

The coherent structure profiles are determined from a 
locally-parallel, linear stability theory. Given the known 
complications that arise when using inviscid theory to com- 
pute eigenvalues and eigenfunctions of the damped por- 
tions of the instability wave, a compressible, viscous ap- 
proach was used for the stability calculations following 
the work of Morris . 12,13 The compressible equations of 
motion are linearized about a locally-parallel flow and the 
compressible parts of the viscous terms in the momentum 
equations are neglected. In the energy equation, both vis- 
cosity and diffusion effects are neglected. After applying 
the modal decomposition for single azimuthal modes, the 
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governing linearized equations become 


dv 

v in 


— 1 

dr 

r r 

d 2 v 

_i_ 

1 dv 

dr 2 

r dr 

d 2 w 

_i_ 

1 dw 

dr 2 

r dr 

d 2 u 

1 du 



dr 2 

r dr 

where A 

2 = a 2 


A + 


A 2 + 


A 2 + 


iau — iMj (lo — aU)p = 0 
n 2 + 1 


(46) 


- *2n -n dp 

v x- w — Re — = 0 

r 2 Or 


n~ + 1 


w - 


z2n „ mRe „ 

— p = 0 


„ -dU „ „ , 

u — p— — n — iaRep = 0 
or 


waves, the linearized equations are the same except the w 
is replaced by — iw in equations (46). 

Using the boundary conditions given by Morris, 13 we 
obtain an eigenvalue problem for a. We used a finite- 
difference approximation to discretize the system of equa- 
tions (46). The eigenvalue was found from diagonalizing 
the resulting matrix using a Newton iteration for refine- 
ment. The four eigenfunctions in these equations were then 
found using the inverse power method and normalizations 
given previously. The density disturbance eigenfunction 
was obtained from the continuity equation. 


—iv 


P = 


(u> - aU) dr 


df> 1 pM 2 p 


(47) 


The eigenvalue depends not only on the flow profiles, 
but it also depends on the parameters n,Mj,u>, and Re. 
The mode number, the Mach number, and the frequency 
are easily chosen parameters of interest. The choice of 
Reynolds number is not as certain. Morris 13 has shown 
how the eigenvalue can depend on the choice of Re. In cal- 
culations of viscous stability, setting Re to a constant has 
been typical. 14 We propose here that since the eigenvalues 
are computed based on local flow conditions that a local 
Reynolds number be used, also. Using the numerical mean 
flow calculations, we determine a local Re from 


PjUjD 
(At) 


max 


(48) 


based on equation (45). This represents the maximum ef- 
fect that turbulent viscosity has on the large structures. 

Though this approach to local stability includes viscous 
effects, it is found to have the same limitations as linear, 
inviscid theory. 15 The analytically continuous eigenvalue 
solution for a damped, supersonically traveling mode has 
corresponding eigenfunctions that diverge for large r rather 
than converging as required to compute the integrals I p and 
I w d- This limits the Mach number for the application of this 
integral technique. 


NUMERICAL RESULTS 

Mean flow 

The calculations were performed for seven jets with exit 
Mach numbers ranging from 0.3 to 2.1 in increments of 



Fig. 1 Jet mean flow profiles computed at 8 axial locations. 
Mj = 1.8; To = T^. 

0.3. In all cases, the jets were cold with the jet exit to- 
tal temperature equal to the ambient temperature. As the 
flow becomes supersonic, the jets were perfectly expanded. 
An example of the computed jet profiles is shown in Fig- 
ure 1 for a Mach 1.8 cold jet. The mean flow profiles have 
smooth gradients in the radial direction and smoothly tran- 
sition from the potential core region to the fully developed 
region downstream. The radial velocity exhibits physically 
realistic behavior. The flow is outward, a positive velocity, 
toward the shear layer from the high speed potential core 
and inward, a negative velocity, outside of the jet indicat- 
ing entrainment of the outer fluid. These flow calculations 
meet the requirements indicated for the stability calcula- 
tions. 

The mean flow calculations also provide the local mo- 
mentum thickness and the local Reynolds number required 
for the stability and the energy integral calculations. The 
momentum thickness, 9, is computed based on only the ve- 
locity profile and the local Reynolds number is computed 
from equation (48) during the mean flow computations as 
part of the turbulence model. The results for the seven 
jets is shown in Figure 2. The jet spreading rate decreases 
with Mach number, as expected, with the momentum thick- 
ness remaining about the same, 9 w 0.22, at the end of 
the potential core as indicated for the Mj = 0.3 and the 
Mj = 1.8 jet cases. The local Reynolds number decreases 
rapidly through the potential core from its initial value at 
the jet exit. 

Local stability 

Using the mean flow as input, the local stability char- 
acteristics were computed for each jet. The eigenvalue or 
local wave number was determined from equations (46) 
where the pair of oblique waves form was used in the com- 
putations. Results are shown in Figure 3 for each jet where 
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Fig. 2 Axial variation of local momentum thickness, 6, and 
local Reynolds number, Eq. (48), computed for 7 jet Mach 
numbers. 

the Strouhal number was set to 0.2 and the azimuthal mode 
number was n = ±1. The local wave number, a, is com- 
plex consisting of a real part, a r , that governs the axial 
phase change of the instability wave and an imaginary part, 
a,:, that relates to the local growth or decay of the wave. In 
this figure, a negative a, represents a growing mode and a 
positive value represents a decaying mode. We can see that 
the growth rate decreases with an increase in Mach number. 
This is expected; but, the decrease in growth rate is also a 
function of the local decrease in Reynolds number. 13 

The local phase velocity is computed from oj/a r . In Fig- 
ure 3, it is shown relative to the ambient speed of sound. 
The relative phase velocity increases directly as the Mach 
number increases. It is an indication of the ability of the 
instability wave to directly radiate noise as the phase ve- 
locity becomes supersonic relative to the ambient speed of 
sound. Low speed jets do not radiate sound from this type 
of source. As the jet speed increases, the higher phase ve- 
locity combined with the growth and decay behavior of the 
instability wave contribute to the ability of these waves to 
radiate noise. 15 

Local energy integrals 

The four integral terms in equation (43) govern the ad- 
vection, transfer, and dissipation of the mean kinetic en- 
ergy. The results of the computations of these integrals for 
the seven jets are shown in Figure 4. As a function of 6, the 
advection of the mean flow kinetic energy, as represented 
by the integral I am , is essentially the same for these jets. 
This indicates that the integral is related to the jet spreading 
rate. The rate of change of I am is initially constant and then 
gradually decreases. The energy moves downstream as the 
jet spreads and then downstream of the potential core, the 
advection slows. 

The next two integrals, I mw \ and I mw 2, represent the 
transfer of energy from the mean flow to the coherent struc- 
tures. The integral I mw i has the classical form for the ve- 
locity shear driven generation of the large-scale structures. 


4.0 

3.5 

3.0 
a 2.5 

2.0 

1.5 
1.0 



0.0 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 

0 



Fig. 3 Local wave number a = a r + iai and phase veloc- 
ity, c p h = cu/a,., as a function of jet momentum thickness, 
n = ±1; fD/Uj = 0.2. 

It is negative indicating energy transfer from the mean flow 
to the modes. Its magnitude is larger for low speed jets and 
decreases in magnitude as the Mach number increases be- 
coming a less efficient energy transfer process. However, 
this view has to be balanced by the presence of I mw 2 . This 
second integral is present due to compressibility effects. At 
low speed, the magnitude of I mw 2 is negligible. Its effect 
steadily increases as the Mach number increases until at 
Mj = 2.1, for example, I mw \ and I mw 2 have similar mag- 
nitude. Part of the reason for the increased magnitude of 
Irnw '2 can be seen by explicitly substituting for p in the in- 
tegral using equation (47). The shear of the density profile 
is a dominant factor in determining the magnitude of I mw 2 ■ 
As the Mach number increased, the jet density increased 
relative to ambient conditions increasing the density shear 
rate. 

One computational difficulty created with the current 
viscous stability method is the presence of the pole in p 
as the imaginary part of a passes from negative to posi- 
tive (Eq. (47)). This must be accounted for in the integrals 
containing p along with the subsequent branch cut associ- 
ated with damped modes. The discontinuities present in the 
plots of I mw 2 are clear evidence of the problem. 

The final integral, I m d, concerns viscous dissipation of 
the mean kinetic energy. In combination with the local 
Reynolds number, it is initially two orders of magnitude 
less than the other integrals and only becomes effective 
when the energy transfer integrals approach zero. The dis- 
sipation is about constant in the potential core region and 
then decreases downstream. There is more dissipation for 
lower speed jets than for the higher speed jets. 

The next set of integrals, shown in Figure 5, govern the 
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Fig. 4 Energy integrals from the mean flow energy equation 
(43). n = ±1; fD/Uj = 0.2. 
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Fig. 5 Energy integrals from the coherent structure energy 
equation (44). n = ±1; fD/Uj = 0.2. 

kinetic energy in the coherent modes; in this case, a pair 
of oblique modes. The energy comes from the mean flow 
through I . mw i which now is positive. As it decreases in 
magnitude, the work done by the coherent mode pressure 
gradients, I p , and the dissipation term begin to reduce the 
energy in the coherent mode. Soon, the dissipation quickly 
dominates to rapidly decrease the mode energy. The effect 
of the compressibility integral I m w 3 has some initial effects 
of adding energy to the coherent modes, but it quickly has 
negligible effects on further energy transfers. 

Nonlinear calculations 

Once the energy integrals were computed, the instability 
wave amplitude function was determined from the solu- 


Fig. 6 Amplitude function magnitude comparisons between 
nonlinear and linear calculations, n = ±1; fD/Uj = 0.2; 

jA 0 j 2 = 10~ 3 . 
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Fig. 7 Shear layer growth comparisons for nonlinear calcu- 
lations. n = ±1; fD/Uj = 0.2; jA 0 | 2 = 10 3 . 

tion of the nonlinear differential equations (43) and (44). 
An initial value differential equation solver implementing 
a degree- 1 Taylor series scheme with fixed point iteration 
was used. 16 The energy integral results were interpolated 
as necessary to meet the requirements of the scheme. Only 
the choice of the initial amplitude was needed to start the 
integration at x = 0, where 6 was given by the flow calcu- 
lations. 

The results for the computed amplitude function mag- 
nitudes for the seven jets are shown in Figure 6. The 
Strouhal number was set to 0.2, the mode number was 
n = ±1, and the initial amplitude |A 0 | 2 was 10 -3 . The 
solid lines represent the nonlinear development of the am- 
plitude function. The dashed lines are the correspond- 
ing linear instability wave amplitude functions determined 
from |A 0 | exp(— f Q a, d\). It can be seen that in this case 
the nonlinear amplitude reaches a peak and begins to de- 
cay much sooner than the linear amplitude. This effect 
becomes more pronounced as the Mach number increases. 
At this level of initial amplitude, the mode energy causes 
increased jet spreading that is more pronounced as Mach 
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Fig. 8 Amplitude function magnitude comparison with and 
without the coherent structure compressibility effects integral 
Imw3 • Mj = 1.8; n = ±1; fD/Uj = 0.2; jA 0 | 2 = 10 3 . 

number increases as shown in Figure 7 for three example 
jet Mach numbers. The increased jet spreading conversely 
lessens energy transfer into the mode leading to earlier 
saturation and decay for higher Mach number nonlinear 
amplitudes compared to linear amplitudes. 

It is interesting to notice that the initial small nonlinear 
amplitude does not coincide with the linear amplitude. In- 
tuitively, the expectation would be that nonlinear behavior 
would begin to take affect after a period of linear growth. 
An examination of the energy exchange integrals is re- 
vealing. As mentioned earlier, the integral has the 

effect of initially adding mean flow energy into the coher- 
ent modes. It contains terms with the radial mean velocity, 
V, and the mean density shear, dp/ dr, which means that 
Imw 3 would not appear in a parallel, incompressible flow 
analysis of this type. For the computed results presented 
here, shown in Figure 5, I m w3 causes the amplitude to 
grow faster immediately and not follow the linear growth. 
To confirm this notion, the calculation used for Figure 6 
was repeated for the Mj = 1.8 jet setting I mw 3 = 0. 
The results are shown in Figure 8. Obviously, without 
the additional energy input into the coherent modes from 
non-parallel flow and compressibility effects, the nonlinear 
amplitude function does initially follow the linear growth 
before nonlinear effects take over. 

The initial amplitude used in the nonlinear calculations 
also affects the results. This is shown in Figure 9 for the 
Mj = 1.8 jet where the amplitude is normalized by its ini- 
tial value. The initial amplitude ranges from |A 0 | 2 = 10 -6 
to 1 0 2 . The normalized linear amplitude line is the same 
for all initial amplitudes. Hence, only one line is necessary 
for comparison. The nonlinear amplitude peaks and be- 
gins to decay sooner as the initial amplitude increases. The 


Fig. 9 Amplitude function magnitude comparisons with dif- 
ferent initial amplitudes. Mj = 1.8; n = ±1; fD/Uj = 0.2. 
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Fig. 10 Shear layer growth comparisons with different initial 
amplitudes. Mj = 1.8; n = ±1; fD/Uj = 0.2. 


stronger initial amplitude leads to the initial increased jet 
spreading (Figure 10) that results in less total energy being 
transferred into the mode from the mean flow. The end re- 
sult is quicker saturation and decay of the mode amplitude. 
This effect was also noted for the two-dimensional, com- 
pressible shear layer. 5 The weaker initial amplitudes result 
in delayed spreading, allowing more energy to transfer into 
the mode even to the point that the amplitude grows larger 
than the linear amplitude. These results point to an issue 
that occurs when comparing computed results to measured 
data. In linear calculations, the unknown initial amplitude 
is used as a parameter to match data with computations. 
With nonlinear calculations, every different initial ampli- 
tude changes the peak location and shape of the amplitude 
function. To illustrate. Figure 11 shows qualitative com- 
parisons of both nonlinear and linear amplitude functions 
to measured data. 17 The data is from a Mach 1.8 cold jet 
and the Strouhal number is 0.25. The measurements in- 
clude a spatial integration and include the presence of all 
mode numbers. The calculations are for mode n = ±1; 
hence, the comparison is inexact. 
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Fig. 11 A comparison of amplitude function magnitude 
calculations versus measured data. 17 Mj = 1.8; n = ±1; 
fD/Uj = 0.25. 

CONCLUDING REMARKS 
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Reynolds -number supersonic jet. J. Fluid Mech., 
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In this method, the large-scale structure is approximated 
by a single -frequency, instability wave of a given azimuthal 
mode number. The fact that the large-structures have a 
spectrum of frequencies and multiple azimuthal modes has 
been ignored. Furthermore, the interactions among these 
modes and with the background, fine-scale turbulence have 
been neglected. In addition, non-parallel flow effects are 
not fully taken into account. 

Despite these approximations, the results are clearly en- 
couraging indicating qualitative agreements with experi- 
ment. The modes initially grow linearly and then satu- 
rate as a result of the nonlinear interactions and the flow- 
divergence effects. Interestingly, the results indicate that 
the initial growth of the modes becomes larger compared 
to that of the linear, parallel-flow case due to non-parallel 
flow and compressibility effects. The results also indi- 
cate that the initial level of the disturbance influences the 
development of the coherent structure and hence the radi- 
ated noise field. As such, experimental data on jet noise 
measurements need to be related to the initial level of the 
disturbances. 
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